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Abstract 

A  two-dimensional  theoretical  model  is  developed  in  this  study  to  simulate  the  transport  phenomena  of  oxygen  in  cathode  gas  diffusion  layer 
(GDL)  of  proton  exchange  membrane  fuel  cell  (PEMFC).  An  analytical  solution  is  then  obtained  accordingly  to  characterize  the  effects  of  GDL 
on  cell  performance.  It  is  found  that  the  concentration  flux  of  oxygen  across  the  GDL  is  primarily  dominated  by  the  thickness  and  porosity  of 
GDL.  For  a  thicker  GDL,  the  diffusion  resistance  increases  and  thus  lowers  the  cell  performance  especially  under  high  current  density  condition. 
On  the  other  hand,  an  increase  of  porosity  will  enhance  the  transport  of  oxygen  and  result  in  significant  improvement  of  cell  performance.  The 
influences  of  system  parameters  including  the  temperature,  channel  height,  inlet  velocity,  and  inlet  pressure  on  the  diffusion  of  oxygen  in  GDL 
are  also  examined  systematically.  Results  provide  insights  into  the  characteristics  of  oxygen  diffusion  in  GDL  and  benefit  the  optimal  design  of 
PEMFC. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

In  the  past  decade,  proton  exchange  membrane  fuel  cell 
(PEMFC)  has  received  much  attention  due  to  its  widely  promis¬ 
ing  applications  in  portable  power  sources,  automobile  power 
systems,  and  stationary  power  plants.  Numerous  studies  have 
devoted  to  develop  theoretical  models  in  order  to  simulate  the 
transport  of  reactants,  electrochemical  reactions,  and  character¬ 
istics  of  membrane  in  PEMFC  [1-22].  Particularly,  since  it  is 
well  known  that  the  performance  of  PEMFC  is  cathode  limited, 
most  of  the  research  efforts  have  focused  on  the  cathode  side 
of  PEMFC  [8-22]  to  investigate  the  composition  and  perfor¬ 
mance  of  catalyst  layer,  transport  of  reactants  in  feed  channel 
and  porous  electrode,  and  water  management  problem. 
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Recently,  many  investigations  [23-3 1  ]  have  concentrated  on 
the  GDL  because  of  its  substantial  influence  on  fuel  cell  per¬ 
formance.  The  GDL  makes  the  gas  reactant  transports  to  the 
catalyst  layer  in  a  diffusive  way  rather  than  direct  contact  in 
order  to  avoid  the  possible  breakdown  of  the  membrane  caused 
by  the  strong  impact  of  reactant  on  the  membrane  surface.  How¬ 
ever,  the  existence  of  GDL  affects  the  mass  transport  of  reactant 
profoundly  and  results  in  a  loss  of  cell  performance.  Thus,  it 
is  quite  important  to  understand  the  characters  of  GDL  for  the 
design  work  of  PEMFC.  Jordan  et  al.  [23]  performed  an  experi¬ 
mental  study  to  find  the  effect  of  carbon  morphology  in  GDL  and 
proposed  a  model  to  explain  how  the  hydrophobicity  and  poros¬ 
ity  of  GDL  affect  water  impregnation  and  gas  diffusion  through 
this  layer.  Park  et  al.  [24]  considered  the  GDL  consists  of  the 
gas  diffusion  medium  (GDM)  and  the  micro-layer,  and  inves¬ 
tigated  the  water  transportation  in  GDM  with  different  PTFE 
contents.  A  novel  design  of  GDL  has  been  developed  by  Chen 
et  al.  [25]  that  insert  a  water  management  layer  between  the 


C.R.  Tsai  et  al.  /  Journal  of  Power  Sources  160  (2006)  50-56 


51 


traditional  GDL  and  the  catalyst  layer.  The  effects  of  different 
fabrication  methods  and  thickness  of  GDL  were  also  discussed 
by  Lee  et  al.  [26].  Mirzazadeh  et  al.  [27]  studied  the  role  of 
GDL  in  the  cell  performance  and  evaluated  the  electrode  with 
or  without  this  layer.  They  found  that  the  effect  of  GDL  is  more 
pronounced  at  higher  current  densities.  Jeng  et  al.  [28  ]  developed 
a  two-dimensional  model  for  oxygen  mass  transfer  in  the  GDL 
of  air-breathing  PEMFC  under  the  influence  of  current  collector 
ribs.  The  resulting  equations  were  discretized  and  solved  by  a 
numerical  technique.  Their  results  indicate  that  the  GDL  is  used 
only  partly  in  the  mass  transfer  process  due  to  the  existence 
of  ribs.  Williams  et  al.  [29]  examined  the  limiting  factors  for 
the  oxygen  transport  in  the  cathode  GDL  with  several  kinds  of 
commercial  GDL.  They  further  studied  the  influence  of  convec¬ 
tion  through  the  GDL  on  the  cathode  limiting  current  by  using 
three  GDLs  with  different  gas  permeability  [30].  A  computa¬ 
tional  analysis  for  the  effect  of  GDM  permeability  on  convective 
transport  also  has  been  performed  by  Pharoah  [31]. 

Although  the  GDL  plays  an  important  role  in  the  optimiza¬ 
tion  of  PEMFC,  it  is  still  difficult  to  evaluate  the  influence  of  an 
individual  property  of  GDL  on  cell  performance  in  experiments 
or  theoretical  analyses,  because  the  electrochemical  reactions 
in  gas  diffusion  electrode  involve  complicated  factors  and  they 
are  strongly  coupled  in  general.  Therefore,  a  simplified  model  is 
necessary  and  helpful  to  understand  the  characteristics  of  GDL. 
In  the  present  study,  we  have  developed  a  two-dimensional  theo¬ 
retical  model  for  the  GDL  along  the  cathode  channel  of  PEMFC. 
The  system  is  simplified  and  thus  an  analytical  solution  can 
be  obtained  to  characterize  the  transport  of  oxygen  in  GDL. 
Accordingly,  the  influences  of  physical  factors  of  the  system 
on  the  cell  performance  can  be  easily  examined  based  on  this 
simple  analytical  solution. 

2.  Theoretical  model  and  formulation 


2.  The  catalyst  layer  is  very  thin  and  can  be  treated  as  a  boundary 
adjacent  to  the  upper  surface  of  GDL.  The  electrochemical 
reaction  occurs  at  this  layer  only. 

3 .  The  GDL  is  isotropic  and  the  transport  of  oxygen  in  this  layer 
can  be  described  by  Darcy’s  law. 

4.  The  generation  and  transport  of  water  are  ignored. 

5.  The  oxygen  in  the  system  behaves  as  an  ideal  gas. 


Accordingly,  the  continuity  equation  in  the  GDL  can  be 
expressed  by 

d(pvx)  d(pvy)  = 

dx  +  dy  ~  {  ’ 


where  p  is  the  density,  vx  and  vy  are  respectively  the  x  and  y 
components  of  the  superficial  velocity  which  is  defined  as  the 
volume  flow  rate  through  an  unit  cross-sectional  area  in  a  porous 
medium.  By  utilizing  the  Darcy’s  law,  the  velocity  components 
vx  and  vy  can  be  written  as 


k  dP  _  K  dP 

/x  dx  ’  y  /x  dy  ’ 


(2) 


where  k  is  the  permeability,  /x  the  dynamic  viscosity,  and  P  is 
the  pressure.  According  to  the  ideal  gas  equation  of  state,  the 
pressure  P  can  be  substituted  by 


P  =  CdRT, 


(3) 


where  Cd  is  the  molar  concentration,  R  the  universal  gas  con¬ 
stant,  and  T  is  the  absolute  temperature.  Note  that  p  =  CdM  in 
which  M  is  the  molecular  weight  of  oxygen.  Eqs.  (1) — (3)  can  be 
combined  together  in  the  form: 


V2 


Cj(x,  y ) 
2 


=  0, 


(4) 


Consider  the  two-dimensional  model  of  cathode  GDL  as 
shown  in  Fig.  1 .  The  GDL  with  thickness  H  is  confined  between 
the  catalyst  layer  and  the  gas  channel.  The  interface  between  the 
GDL  and  the  gas  channel  is  located  at  y  =  0.  Pure  oxygen  is  fed 
into  the  gas  channel  with  inlet  velocity  vq,  where  the  channel 
height  and  length  are  hc  and  L,  respectively.  To  study  the  charac¬ 
teristics  of  transport  and  consumption  of  oxygen  in  this  system, 
some  simplifications  are  made  as  list  below: 

1 .  The  system  is  isothermal  and  in  a  steady  state. 


_ -*  Catalyst  layer 

Las  diffusion  layer  (GDL) 

"  v0 

. 

Gas  channel 

hc 

x= 

Bipolar  plate 

Fig.  1.  Schematic  representation  of  the  two-dimensional  model  for  the  gas  dif¬ 
fusion  layer  in  cathode  side  of  PEMFC. 


where  V2  =  d2/dx2  +  d2/dy2  is  the  Laplace  operator.  Define: 


0(x,  y)  = 


Cj(x,  y) 

2 


(5) 


thus,  Eq.  (4)  becomes 

V2[0(.r,  >■)]  =  0.  (6) 


The  above  equation  characterizes  the  concentration  distribution 
in  the  GDL  and  can  be  solved  simply  with  appropriate  boundary 
conditions.  At  the  side  walls  of  GDL  x  =  0  and  x  =  L,  there  is  no 
concentration  flux  and  the  boundary  conditions  are 


d& 

dx 


x=0 


80 

dx 


=  0. 

x=L 


(7) 


At  the  interface  between  the  GDL  and  the  gas  channel  y  =  0,  the 
concentration  is  continuous  across  the  interface  and  the  bound¬ 
ary  condition  is  given  by 


C1 

0 1  y=o  =  exp  |  -  — 


2x  \ 

~L~J  ' 


(8) 


52 


C.R.  Tsai  et  al.  /  Journal  of  Power  Sources  160  (2006)  50-56 


In  above  equation,  the  results  for  the  profile  of  oxygen  concen¬ 
tration  along  the  channel  given  by  Dohle  et  al.  [11]  are  used  and 
the  corresponding  equation  is 

Cd(*,0)  =  C0exp^--^  ,  (9) 

where  Co  is  the  inlet  concentration  and  L\  the  characteristic 
length  of  oxygen  consumption  at  the  reaction  order  y  =  1 .  To 
derive  the  boundary  condition  at  the  interface  between  the  GDL 
and  the  catalyst  layer  y  =  H,  we  can  define  the  mass  flux  of  oxy¬ 
gen  wo  as  a  function  of  local  current  density  i(x)  according  to 
the  electrochemical  reaction  at  the  catalyst  layer: 

02  +  4H+  +  4e“  2H20,  (10) 

so  we  have: 


and  then  Eq.  (15)  can  be  solved  to  get  the  boundary  condition  at 
y  =  H: 


C2(x)G2 
~  2 
C2(x) 

=  — exp 


/  2(e  —  1)  ioH 

V  neDlF  Cref  6XP 


(17) 


where 


G  —  exp 


/(e-  1)  i0H 
\HSDqF  Cref  6XP 


(18) 


Eq.  (6)  together  with  the  boundary  conditions  (7),  (8),  and  (17) 
can  be  solved  to  obtain  an  analytical  solution  in  the  form: 


&(x,  y)  = 


1  —  exp 


y 


mo\y=H  =  PVy  =  MCdVy  = 


iM 

nF’ 


(11) 


where  n  =  4  is  the  number  of  electrons  transferred  per  mole  of 
oxygen  consumed  and  F  the  Faraday  constant.  The  mass  flux 
m  o  also  satisfies  the  relation  [9]: 


n=  1 


1  —  exp  ( —  2(  j  cos(«7r) 


kl  (TL\2  4-  JL 

2  \  L  )  'T  Li 


X  cos 


fmt  \  j  G2  sinh^y 
\~LX)  |  sinh 


mo\y=H  =  -£DeQM— - 

dy 


y=H 


FsMvyCd  | 


(12) 


[sinh  ('f  y)  -  tanh  (f  H)  cosh  (f  y)] 
tanh  (fH) 


.  (19) 


where  the  first  and  second  terms  on  the  right  hand  side  are 
the  amounts  of  oxygen  flux  caused  by  diffusion  and  convec¬ 
tive  motion,  respectively.  Here,  s  is  the  porosity  of  GDL  and  D q 
the  effective  molecular  diffusion  coefficient  of  oxygen  that  is 
related  to  the  diffusion  coefficient  in  a  nonporous  domain  Dq  by 
Dq  =  Doe1'5  [9].  By  equating  both  Eqs.  (11)  and  (12),  we  may 
obtain: 


3Cd  £  —  1 

-  =  - i\x) 

9y  y=H  hsFDq 


(13) 


The  Tafel  equation  [9,12]  is  used  to  describe  the  distribution  of 
current  density  along  the  boundary  of  catalyst  layer: 


i(x)  —  i  o 


/  Cd|y=ff\ 

V  cref  J 


(14) 


where  io  is  the  exchange  current  density,  Cref  the  reference  con¬ 
centration  of  oxygen,  a  the  charge  transfer  coefficient  and  rj  the 
electrode  activation  overpotential.  Therefore,  Eq.  (13)  becomes 


dCd 


9v 


y—H 


(e  -  1)  /o  ( uF  \ 
HSFDq  Cref  EXP  \RT  11 ) 


Cd\y=H- 


(15) 


Note  that  Eq.  (9)  indicates  the  concentration  of  oxygen  decays 
exponentially  along  the  channel.  Hence,  it  is  reasonable  to 
assume  that  the  concentration  of  oxygen  in  the  GDL  also  behaves 
in  a  similar  way.  Thus,  the  concentration  function  Cd(x,  y)  can 
be  written  by  separation  of  variables: 


(16) 


This  equation  gives  the  concentration  of  oxygen  in  the  GDL  in 
a  simple  analytical  form. 

To  determine  the  voltage-current  curves  of  the  PEMFC,  we 
define  the  average  current  density  i  over  the  channel  length  L  as 

1  /'L- 

i  =  I  i(x)  dx,  (20) 

L  J  o 

and  the  operating  cell  voltage  is  given  by 

17  =  Vcell  ri  *)ohmi  (21) 


where  Vceii  is  the  open  circuit  voltage,  ;;  can  be  determined 
from  the  Tafel  Eq.  (14),  and  ^()|lm  is  the  ohmic  losses.  Since 
in  most  PEMFCs  the  electrical  resistance  is  mainly  caused  by 
the  ion  transport  through  the  membrane,  the  ohmic  losses  are 
approximated  by  a  plain  expression: 


^ohm  — 


;7m 

°m 


(22) 


where  /m  is  the  membrane  thickness  and  erm  the  membrane  con¬ 
ductivity.  Eqs.  (19)-(22)  provide  a  simple  technique  to  estimate 
the  influence  of  GDL  on  fuel  cell  performance. 


3.  Results  and  discussion 


We  first  examine  the  influence  on  cell  performance  due  to 
the  existence  of  GDL  and  compare  the  results  with  those  of 
Dohle  et  al.  [11]  in  which  they  considered  the  case  without 
GDL.  Fig.  2  illustrates  the  variations  of  I-V  curves  for  sev¬ 
eral  assigned  values  of  the  GDL  thickness.  Based  on  the  same 
parameters  used  in  Ref.  [11],  it  is  found  that  the  present  results 


Cd(x,  y)  =  C(x)9(y), 
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Current  density  (A/cm2) 

Fig.  2.  Voltage-current  curves  for  four  assigned  values  of  GDL  thickness. 


are  in  agreement  with  those  of  Ref.  [11]  for  the  case  without 
GDL.  The  use  of  GDL  apparently  lowers  the  cell  performance 
and  the  curve  falls  gradually  with  increase  of  GDL  thickness. 
This  effect  is  more  pronounced  especially  at  high  current  den¬ 
sity  condition,  which  indicates  that  a  thicker  GDL  gives  rise  to 
a  higher  resistance  for  oxidant  transport  toward  the  electrode 
surface.  The  influence  of  GDL  for  different  channel  length  is 
demonstrated  in  Fig.  3  for  three  assigned  values  of  L.  It  has  been 
proposed  by  Dohle  et  al.  [11]  that  a  decrease  of  channel  length 
will  improve  cell  performance  in  the  whole  range  of  current 
densities  because  the  oxygen  concentration  is  more  uniform  in 
shorter  L  and  this  effect  raises  the  overall  performance  per  unit 
area  as  shown  in  Fig.  3  for  the  cases  without  GDL.  When  we  take 
the  GDL  into  consideration,  it  is  found  that  the  influence  of  GDL 
is  quite  significant  especially  in  short  channel  length.  For  the 
case  L  =  400  cm,  the  limiting  current  density  reduces  only  11% 
approximately.  While  for  the  case  L=  100  cm  it  reduces  about 
29%,  which  reveals  that  it  is  very  important  to  improve  the  losses 
caused  by  the  use  of  GDL  particularly  in  relatively  short  channel 
length. 

Table  1  gives  the  values  of  parameters  used  in  the  follow¬ 
ing  analysis  for  a  base  operating  condition  of  the  fuel  cell. 


Fig.  3.  Volage-current  curves  for  three  assigned  values  of  channel  length  L  with 
and  without  GDL. 


Table  1 


Values  of  physical  parameters  used  in  the  base  case 


Cell  temperature,  T  (K) 

343 

Inlet  pressure,  P  (atm) 

1 

Inlet  flow  velocity,  vo  (cms-1) 

145 

Inlet  concentration,  Co  (mol  cm-3) 

3.87  x  1(T5 

Reference  concentration,  Cref  (mol  cm-3) 

3.18  x  10~5 

Channel  length,  L  (cm) 

44 

Channel  height,  hc  (cm) 

0.1 

Membrane  thickness,  /m  (cm) 

0.018 

Gas  diffusion  layer  thickness,  H  (cm) 

0.04 

Transfer  coefficient,  a 

1 

Number  of  electrons,  n 

4 

Exchange  current  density,  *o  (A  cm-2) 

5.74  x  10~7 

Molecular  weight,  M  (kg  mol-1) 

32 

Dynamic  viscosity,  /x  (gem-1  s-1) 

2.03  x  10~4 

Permeability,  k  (cm2) 

10~8 

Universal  gas  constant,  R  (J  mol-1  K-1) 

8.314 

Faraday  constant,  F  (C  mol-1) 

96487 

Porosity,  s 

0.4 

Diffusion  coefficient  of  oxygen.  Do  (cm2  s-1) 

0.0991 

Ideal  open  circuit  voltage,  Vceii  (V) 

1.2 

Conductivity,  am  (1  £2-1  cm-1) 

0.0868 

The  variations  of  the  distribution  of  oxygen  concentration  in 
the  gas  channel  and  GDL  for  three  assigned  values  of  poros¬ 
ity  s  under  the  same  output  of  cell  voltage  0.5  V  are  shown  in 
Fig.  4(a)-(c)  and  the  corresponding  I-V  curves  are  shown  in 
Fig.  5. The  channel  height  hc  and  the  thickness  of  GDL  H  are 
respectively  0. 1  and  0.04  cm  with  the  interface  at  Y=  0. 1  cm 
indicated  by  the  dashed  line.  One  can  see  that  at  e  =  0.3  in 


Fig.  4.  Distribution  of  oxygen  concentration  in  the  cathode  gas  channel  and 
GDL  for  three  assigned  values  of  porosity  with  the  same  output  of  cell  voltage 
0.5  V;  (a)  e  =  0.3,  (b)  e  =  0.4,  (c)  £  =  0.5. 
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Fig.  4(a),  the  oxygen  concentration  within  the  GDL  descends 
rapidly  toward  the  electrode  surface  and  most  of  oxygen  near 
the  catalyst  layer  is  consumed  due  to  the  low  diffusivity  in 
accord  with  the  low  porosity  of  GDL.  Hence,  the  concentra¬ 
tion  flux  through  the  GDL  is  small  and  fails  to  supply  sufficient 
oxidant  to  the  catalyst  layer  for  electrochemical  reaction.  Conse¬ 
quently,  the  cell  voltage  falls  quickly  against  the  current  density 
as  shown  in  Fig.  5.  As  the  porosity  of  GDL  increases,  for  exam¬ 
ple,  e  =  0.4  and  0.5  as  shown  in  Fig.  4(b)  and  (c),  the  mass 
transport  through  the  GDL  will  be  increased  and  therefore  the 
voltage  loss  due  to  the  presence  of  GDL  will  diminish  obviously 
as  indicated  in  Fig.  5.  These  results  suggest  that  the  porosity  of 
GDL  plays  an  important  role  in  the  factors  to  affect  the  cell  per¬ 
formance  and  should  be  as  large  as  possible  in  practical  design  of 
PEMFC. 

The  operating  temperature  of  fuel  cell  is  also  an  important 
factor  and  its  effect  on  the  mass  transport  of  GDL  is  displayed  in 
Fig.  6(a)-(c)  via  the  velocity  profiles  in  the  GDL  for  three  given 
temperatures.  It  is  noted  that  actually  many  parameters  in  the 
system  are  functions  of  temperature  and  the  influences  of  tem¬ 
perature  variation  on  cell  performance  are  quite  complicated. 
Here,  we  simply  consider  the  effect  of  temperature  arising  from 
the  ideal  gas  equation  as  indicated  in  Eq.  (3)  and  assume  the 
other  parameters  are  invariant  with  temperature  to  simplify  the 
analysis.  At  lower  operating  temperature  as  shown  in  Fig.  6(a) 
for  T  =  323  K,  it  is  found  that  the  x-component  of  velocity  for 
the  gas  in  GDL  is  very  small  and  the  direction  of  velocity  is 
almost  toward  and  perpendicular  to  the  catalyst  layer.  Thus,  the 
diffusive  direction  of  oxygen  in  GDL  is  primary  dominated  by 
the  catalyst  layer  where  the  electrochemical  reaction  occurs. 
Simultaneously,  the  magnitude  of  velocity  decays  along  x-axis 
direction  while  grows  gradually  approaching  the  catalyst  layer. 
Similar  results  can  also  be  observed  in  Fig.  6(b)  and  (c)  for 
T  =  343  and  363  K,  respectively.  Apparently,  an  increase  of  tem¬ 
perature  enhances  the  diffusive  velocity  in  the  GDL  and  thus 
improves  the  cell  performance  as  shown  in  Fig.  7.  While  the 
concentration  is  in  inverse  proportional  to  the  temperature  for 
a  constant  inlet  pressure,  the  improvement  of  cell  performance 
by  increasing  temperature  is  limited.  Moreover,  the  I-V  curve 
for  higher  temperature  appears  to  possess  smaller  limiting  cur- 


Fig.  6.  Velocity  profiles  in  the  GDL  for  three  assigned  temperatures;  (a) 
T=  323  K,  (b)  T=  343  K,  (c)  T=  363  K. 

rent  density  because  the  potential  loss  due  to  mass  transport 
is  more  significant  at  higher  temperature  that  causes  the  corre¬ 
sponding  I-V  curve  fall  faster  early  than  the  lower  temperature 
case. 

The  geometric  effect  of  channel  size  is  a  critical  factor  as  well. 
Here  we  simply  consider  the  variations  of  the  channel  height  and 
the  results  are  demonstrated  in  Fig.  8  for  three  assigned  values 
of  channel  height  with  the  same  inlet  velocity  200cms_1.  It 
is  found  that  the  limiting  current  density  will  be  enlarged  by  an 
increase  of  channel  height  and  the  cell  performance  will  be  raised 
significantly  at  high  current  density.  Since  the  mass  flow  rate  is 
larger  for  a  higher  channel  height,  the  average  concentration  will 
increase  within  the  same  channel  length  and  results  in  better 
cell  performance  especially  at  high  current  density  condition. 


Fig.  7.  Variation  of  voltage-current  curves  with  temperature. 
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Currentdensity  (A/cirf ) 

Fig.  8.  Variation  of  voltage-current  curves  with  channel  height. 


Current  density  (A/cm") 

Fig.  9.  Variation  of  voltage-current  curves  with  inlet  pressure. 


Similar  results  also  can  be  obtained  for  different  inlet  velocities 
with  a  fixed  channel  height  due  to  the  same  causes. 

Fig.  9  illustrates  the  variations  of  I—V  curves  with  three 
different  inlet  pressures.  From  Eq.  (3)  it  is  known  that  the 
concentration  is  proportion  to  the  pressure  for  an  isothermal 
system.  Consequently,  the  cell  performance  appears  to  be  quite 
sensitive  to  the  inlet  pressure  and  limiting  current  density  is 
also  approximately  proportion  to  the  inlet  pressure.  This  find¬ 
ing  suggests  that  it  should  be  an  efficient  way  to  enhance  the 
cell  performance  by  utilizing  a  high  inlet  pressure  in  compar¬ 
ison  with  the  other  factors  such  as  temperature  and  channel 
height. 

4.  Conclusions 

We  have  developed  a  simple  two-dimensional  model  for 
the  mass  transport  of  oxygen  within  the  GDL  of  PEMFC  and 
accordingly  obtained  an  analytical  solution  for  the  distribution  of 
oxygen  concentration  in  GDL.  Based  on  this  analytical  solution, 
the  influence  of  GDL  on  cell  performance  has  been  examined 
and  the  results  indicate  that  the  utilization  of  GDL  indeed  pro¬ 
duce  significant  potential  losses.  Such  losses  can  be  improved  if 
one  used  a  thinner  GDL  or  increases  the  porosity  of  GDL,  and 


the  effect  of  porosity  seems  to  be  more  pronounced.  The  other 
physical  factors  including  the  temperature,  channel  height,  inlet 
velocity,  and  inlet  pressure  are  also  investigated.  It  is  found  that 
the  most  efficient  way  to  enhance  the  cell  performance  is  to 
increase  the  inlet  pressure  since  the  concentration  flux  through 
the  GDL  to  the  catalyst  layer  will  greatly  increase  especially 
under  high  current  density  condition. 

Although  the  system  is  simplified  by  some  critical  assump¬ 
tions  such  as  to  treat  the  catalyst  layer  as  a  thin  boundary 
and  neglect  the  water  transport  problem,  the  present  model 
still  provides  a  simple  and  direct  technique  to  estimate  the 
influence  of  GDL,  and  will  be  helpful  for  the  design  work  of 
PEMFC. 
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